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Beam asymmetries result in statistically anisotropic cosmic microwave background (CMB) maps. 
Typically, they are studied for their effects on the CMB power spectrum, however they more closely 
mimic anisotropic effects such as gravitational lensing and primordial power asymmetry. We discuss 
tools for studying the effects of beam asymmetry on general quadratic estimators of anisotropy, 
analytically for full-sky observations as well as in the analysis of realistic data. We demonstrate this 
methodology in application to a recently detected 9<r quadrupolar modulation effect in the WMAP 
data, showing that beams provide a complete and sufficient explanation for the anomaly. 



The cosmic microwave background (CMB) is a power- 
ful probe for both modern and future cosmology: its rota- 
tionally invariant power spectra have been instrumental 
in hammering out the details of the current concordance 
model, its non-Gaussianities have the potential to dis- 
criminate between various early-universe models, and its 
statistical anisotropies can be used to probe astrophys- 
ically interesting secondary effects such as gravitational 
lensing. 

We observe the CMB through the convolution of an 
instrumental beam, an effect which must be carefully 
treated in analysis. Qualitatively, the effects of beams 
are twofold: (i) they suppress structures on scales smaller 
than the beam size; and (ii) if the beams are asym- 
metric, they can introduce statistical anisotropies into 
the observed CMB which can bias estimators for other 
anisotropic signals. The purpose of this paper is to col- 
lect results on the simulation of beam effects, and to 
present fast, accurate techniques for forecasting and cor- 
recting the effects of beams on estimators of statistical 
anisotropy. In Sec. [I] we present a model of beam asym- 
metries and we discuss the covariance which beams pro- 
duce in the observed CMB in SecHH In Sec. lHII we derive 
the effects of the anisotropic covariance on anisotropy es- 
timators, and in Sec. |IV| we illustrate this discussion by 
applying these techniques to study the effects of beams 
on the highly significant quadrupolar modulation effect 
in the WMAP data. Our conclusions are collected in 
Sec.|V| The effect of beam asymmetries on the estimated 
power spectrum of the CMB for general survey geome- 
tries is discussed in Appendix |B"} 



I. MODEL 

In a realistic CMB observation, the effective sky signal 
at each point in the time-ordered data (TOD) is a convo- 
lution of the true sky signal with the experimental beam, 
oriented according to the scan strategy. Schematically, 



we have 

T t = / d0ri(0)6(0) + rii (1) 
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where Tj is the temperature for time-step i in the TOD, 
0(f2) is the underlying CMB signal, ri(f2) is the beam 
response and n.; is the instrumental noise. For the pur- 
poses of compact notation we will abbreviate J s2 dQ as 
J for the remainder of this paper. The integral in Eq. [l] 
can be performed by brute force in real-space using inter- 
polation on pixelized maps of the beam and sky [T] . For 
this approach to be computationally feasible, the beam 
must be assumed zero outside some small patch where its 
response is peaked, and so it is difficult to study sidelobe 
effects with this approach, although it can be quite fast. 

In this work, we will find it more useful to work in 
harmonic space, where the effects of beams are easier to 
study analytically. We begin by writing the beam re- 
sponse as a harmonic sum. If we center the beam at the 
north pole, with some fiducial beam axis aligned along 
the +x-axis (the <f> — meridian) , and expand it in spher- 
ical harmonics bi m , then the beam at location f2 for the 
jth TOD observation is given by (e.g. [5]) 

S max 'max I 

n(n)= £ ^ D ™{<t>iA,ai)b ls0 Y lm {n). 
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(2) 

For the purposes of compact notation we will drop the 
summation limits in what follows. The limits themselves 
will be discussed later. The action of the Wigner-D ma- 
trix can be visualized as follows: imagine fixing the co- 
ordinate system in space and performing right-handed 
rotations of the beam image about the z-axis by an an- 
gle on, then about the y-axis by an angle Qi, and finally 
about the z-axis by an angle <f>i. The first rotation gives 
the beam its orientation: aj is the angle of the fidu- 
cial beam axis, measured from the southern side of the 
meridian which passes through the pixel location (9i, (f>i) 
assigned to the observation. We use s Yi, n to denote a 
spin- weighted spherical harmonic, of which the standard 
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spherical harmonics are a special case with s = 0. Un- 
less otherwise noted, the harmonics should be taken as 
functions of fi. We can then rewrite Eq. ([I]) as 
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In the second step we have used the close relationship be- 
tween Wigner-Z) functions and the spin-weighted spheri- 
cal harmonics [3] , and introduced the beam transfer func- 
tion Bi s given by 
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For a beam which is normalized to have unit response to 
a monopole, -Boo = 1- Wc shall refer to the s = coef- 
ficients of the beam as the symmetric part, as they rep- 
resent the component of the beam which depends only 
on the radial distance from its center. The s^O coef- 
ficients encapsulate beam asymmetry. On scales much 
smaller than the beam size, Bi s becomes very small, 
which effectively band- limits Bi s ®i m . This scale deter- 
mines Zmax- The evaluation of Eq. ([3| may then be per- 
formed in 0(s max ^ lax ). To perform this convolution for 
each timestep in the TOD is in general prohibitively ex- 
pensive, and some approximations must be made. There 
are several possible approaches: 

(1) The convolution can be computed over a grid cover- 
ing the full rotation group with fast-Fourier-transforms 
for the Euler angles <f> and a and, optionally, for 9 0]. 
The TOD can be obtained by interpolation off this grid. 

(2) If the beam may be represented using a small num- 
ber of symmetric basis functions, each of these may be 
rapidly convolved in harmonic space and then sampled 
based on the location and orientation of these basis func- 
tions for each sample in the TOD [5] . The difficulty here 
is the ability to represent the beam as a sum of symmet- 
ric functions. Note that in the limit that the beam is 
represented as a sum of delta functions, this approach is 
conceptually the same as real-space integration. 

We note the above approaches for completeness. In 
this work, we will use a popular map-based ap- 

proach based on the assumption that the TOD noise is 
uncorrelated on the timescales which separate pixel vis- 
its. In this case, it is a good approximation to the map- 
making process (in the absence of beam deconvolution 
[10] ) to take 



e(fi p ) + n(a 



(5) 



where O is an effectively observed sky, n(Q p ) is a noise 
map, and H p is the number of elements in the sum, which 
is taken over all hits assigned to pixel p, with center 
at dp. This approach can also be used for differencing 
experiments, which suppress correlations between TOD 



samples by mapmaking from the difference between two 
nearly-identical detectors, to remove common mode fluc- 
tuations. In this case, one can use an effective beam 
which is a hit-weighted sum of the two beams which are 
differenced [TT]. 

In conjunction with Eq. the sum of Eq. Q can be 
seen to effect a Fourier transform of the distribution of 
orientation angles, with an effective observed sky given 
by © 
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where the details of the scan strategy are contained in 
the spin — s field 



w(Q p ,-s) 
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Since a Y\ m involves s (spin-weighted) derivatives of the 
Yi m , each term in the s sum is the real-space product 
of the scan strategy and beam-filtered derivatives of the 
CMB. For a beam which is approximately azimuthally 
symmetric, or a scan strategy which broadly distributes 
the orientation angles, w(fl p , —s)Bi s falls off sharply with 
s, and it follows that calculation of only the lowest s 
terms suffice to give a good approximation to the beam- 
convolved map. This determines an effective s max which 
can be much less than that naively required to describe 
accurately the beam in Eq. 

Given a scan strategy, Eq. Q provides an 0(s max Z^ ax ) 
method to compute effectively the sky observed by an ex- 
periment with an asymmetric beam Bi s and given scan 
strategy w(il p , s). This approximation is useful not only 
for its speed, but also to gain an intuitive analytical un- 
derstanding of beam effects, which we proceed to discuss 
in the following sections. 



II. COVARIANCE 

Beam effects are linear in the underlying CMB, and 
so do not affect its (assumed) Gaussianity. For Gaus- 
sian models, the statistics of the observed CMB remain 
completely characterized by its covariance. The effect 
of beams is simply to introduce statistical anisotropics 
which give off-diagonal and m-dependent contributions 
to the covariance. 

In harmonic space the beam-convolved sky is given by 

©i'ro' = ^2 ^ Bis ©im -S W LM 
LMS Im 



x / -sYlaiqY^, sYi m , 



where 



swlm = I w(fl, S) S Y£ M 
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are the spin-,? multipoles of w(tt, S). 

The covariance of the beam-convolved CMB may then 
be written as 
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The A; imi / 2m2 term contains the part of the covariance 
which is due to beam asymmetries. We further split it 
into two terms, such that A = A^ + A^ 2 ). The AW 
terms are those which couple an s = mode of the convo- 
lution with an s 7^ mode. The notation arises because 
we think of them as being first order in any beam asym- 
metry: 

A Lm 1 l 2 m 2 =J2J2 sWLM 
S^O LAI 



x / S Y, 



S*LM 



)Y hmi sYi 2m2 Bi 2S B h oC® e l(lo2) 



(11) 



The (1 -H> 2) represents the interchange of lx>h an d 
mi, m,2 in the preceding expression. The A^ 2 ) terms cou- 
ple two s^O modes: 



Si^O Im 
S 2 ^0 
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Traditionally, CMB analyses have focused on the 
power spectrum, or average diagonal elements of the co- 
variance matrix for each I. The A^ term evaluates 
to zero for these elements, however, and so the effects 
of beam asymmetries on the power spectrum are due 
solely to the terms. For estimators of statistical 

anisotropy, however, the dominant contributions to A 
are expected to be the AW terms, which involve only 
one power of -Bz( s ^o)- Beam asymmetries therefore have 
quantitatively different effects on the power spectrum 
and on anisotropy estimators. A further important dif- 
ference is that beam asymmetries are only an issue for 
the power spectrum on small scales, but, due to mode 
coupling, they can still be important when reconstruct- 
ing large-scale modes of any statistical anisotropy. Note 
that the above arguments only apply exactly on the full- 
sky with uniform pixel weighting. For pseudo-Q power 
spectra on a cut sky or with anisotropic weighting (for ex- 
ample, to mitigate inhomogeneous pixel noise), the beam 
anisotropies can couple with the asymmetry introduced 
by the weights, which gives A^ a contribution to the 
power spectrum. However, this is only significant near 



strong inhomogeneities in the pixel weights, and is there- 
fore generally suppressed. For the remainder of this pa- 
per, we will focus on the effects of beams on anisotropy 
estimators. The effects on the power spectrum are dis- 
cussed in Appendix IB] see also [TT] . 



III. ANISOTROPY ESTIMATORS 

The CMB is assumed to be statistically isotropic to a 
good approximation, but there may be small contribu- 
tions to the covariance from a variety of effects, such as 
gravitational lensing (if the lensing potential is consid- 
ered as fixed, see e.g. [12] for a review), inhomogeneous 
reionization (if the reionization history is fixed, see e.g. 
|13j). Doppler modulation [14], or more exotic statistical 
anisotropy (e.g. |15j and references therein). Following 
the notation introduced in the previous section, we will 
write the CMB covariance as 



hmi l 2 ~m 2 ' 



.(13) 

where i labels the various physical effects that contribute 
to the anisotropy. 

If we assume that the anisotropy from each effect i is 
sourced linearly by multipoles s^Ilm w ith spin- weights 



{&} that satisfy [<?*/£!]* = (-l) s '+ M _ s; £ u , then 
covariance under rotations (i.e. the requirement that 
if is rotated, Jlai must rotate in tandem) and par- 
ity (s'/[']f -> (-l) S ' +i -s>/Sf) generally requires that 
each term in Eq. |l3|) has the form 
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(14) 

where W« ai (I Xj Z 2) L) = (-l) S 'W% _ si (hJ 2 ,L) is a 
weight function which describes the way in which the 
anisotropy field couples the multipoles, while si and 
S2 = S % — si label different partitions of the spin 
S l between two spin- weighted harmonics . Typically, 



1 Two examples may help to solidify the notation. For lensing, 
— > + d a V a to first order in the lensing deflection d a . This 
gives a covariance with nonzero weights 



h 



h 

2^2^2 + 1)' 



and {sfij3} = {±i^LAf } are the spin-weight multipoles of the 
deflection field. Note that we have not assumed that d a is a 
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one of si or s 2 is zero since we are dealing here with 
the spin-0 temperature. It can be seen that the A^ 1 ) 
term of the beam covariance in the previous section is 
of this form (although the smaller A^ 2 ) term is not). 

The Al^^, i are symmetric under (1 «-» 2) so wc 



hmi hm 2 



may take WgJ (h, h, L) = W^} (h,h, L )- Moreover, 
(_1)"»2 aW }a Hermitian under (1 «-> 2) which 

gives rise to the spin-flip symmetry 

W^ si (h,l 2 ,L) = (-lfw%_Jl 2 ,h,L). (15) 

Note that the quantity W{1\, l 2i L)s/lm is essentially 
equivalent to the bipolar spherical harmonic coefficients, 
Aff 1 , of Hajian and Souradeep [IB]. The formalism 
which we will use here and the bipolar spherical har- 
monic formalism can be thought of as two different rep- 
resentations of the same symmetry relations, analogous 
to e.g. Clebsch-Gordan coefficients and Wigner 3j sym- 
bols. Similarly, their relative benefits depend on the use. 
For calculations, the quadratic estimator approach often 
results in simpler expressions, but the A^J^ can prove 
more useful for blind searches, or for gaining insight into 
the physical interpretation of the anisotropies and the 
relationships between different models [TT] . 

As discussed in |15j . optimal quadratic maximum- 
likelihood estimators can be constructed for sJlm, un- 
der the assumption that their effects are perturbative. 
The estimators approximately maximize the Gaussian 
log-likelihood C with respect to the s/lm, so that they 
solve dC/dsfi,M = 0. The solution is constructed from a 

(i) 

set of quadratic building blocks S i h LM , each of the form 
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where 0; m is the inverse-variance filtered observed sky 
(in general the observed sky premultiplied by the signal- 
plus-noise inverse covariance including all anisotropic 
contributions to the covariance that do not depend on 
the set of parameters that are being estimated). The 
inverse-variance filtering can be performed quickly using 
conjugate descent with a good preconditioner, the best 
to date being that of [BJ . The estimator for gi f^ M is then 
given by 



where the ensemble average is taken over realizations of 
the CMB and noise. The "mean-field" term (s^lm) sub- 
tracts off anisotropy due to anisotropic noise, sky cuts, 
and known anisotropic components of the covariance (e.g. 
beam asymmetries). The matrix J 7-1 is the inverse of the 
Fisher matrix, which is calculated as 
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gradient here. For beam asymmetries with sI^lm™*^ 
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where C 1 is the inverse covariance matrix used to con- 
struct Qim- The Fisher matrix can be shown to equal 

(i) 

the covariance of the S ih LM : 



J 7 iS'LMJSiL'M' — (s^LMSJ ^L'M') ~ (s^LAf) (s^i' M') ■ 

(19) 

For an observation and underlying CMB that are statisti- 
cally isotropic (e.g. full-sky coverage with homogeneous 
noise levels and symmetric beams), rotational symme- 
try requires that the inverse-variance filtered CMB has a 
diagonal covariance 

(6^^) = ^ (ko.), (20) 

h 

where l/Cf ot is the inverse- variance filter. This propa- 
gates to the Fisher matrix J 7 , which is then also diagonal 
in L, and independent of M: 
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(21) 



The isotropic Fisher matrix is chiefly useful for forecast- 
ing purposes. In practice, inhomogeneous sky coverage 
and foreground cuts mean that it should be estimated 
from simulations. 

If the weights are separable, such that W{l\,l 2 ,L) = 
Wx(Ii)W 2 {1 2 )Wl(L), or can be decomposed as a sum of 
separable terms, then these estimators have fast position- 
space forms with computational cost C(^, 3 nax ), and the 
isotropic Fisher matrix can be evaluated in O(Zmax) |13) . 

It can be seen clearly from Eq. (Ill that these tools 



apply to the covariance produced by beam asymmetries. 
Optimal estimators could, for example, be formed to re- 
construct the components of swlm for each S. In prac- 
tice, the instrumental scan strategy is fixed and asym- 
metric beams act as a source of bias for other anisotropy 
estimators, which have the form of Eq. ( 16 ). In this view, 



beams simply make a contribution to the covariance of 
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the observed sky. They should be incorporated into the 
inverse-variance filtering operation and the mean-field 
subtraction. For realistic observations, this can be done 
easily as the inverse-variance filtering step is done using 
conjugate descent, and requires only a fast method to 
apply the beam effects, such as that provided by Eq. ([8]). 
This is demonstrated at the TOD level in jTD] , for exam- 
ple. The mean-field can be determined straightforwardly 
from simulations. Analytic calculations are also feasible 
if the inverse- variance filter is isotropic, and can be useful 
for forecasting purposes. Neglecting a known source of 
anisotropy in the data during mean-field subtraction will 
generally bias anisotropy estimators for other effects. Ex- 
plicitly, the mean-field bias on an estimator S j f?L with 
weight function Wg) (h, l 2 , L) by a contaminant i with 
covariance as in Eq. ( |M| is 

\S3 n LM/ - S'JLM ^ I 87r £»tot£<tot 

l 1 L 2 S's 1 s' 1 ( i '2 

x<;,(M 2 ,i) (_';,, 41). (22) 

A nice feature of this result is that the bias rather di- 
rectly traces the spatial distribution of the contaminant: 
the bias is an isotropically filtered version of the con- 
taminant. In the case of pixel- uncorrelated anisotropic 
instrumental noise, for example, the mean-field simply 
traces the anisotropic part of the noise variance map [TH] • 
In the case of beams, on the other hand, the mean-field 
traces the components of the scan strategy swlm- 



IV. CASE STUDY: PRIMORDIAL POWER 
ASYMMETRY 

As an example of this machinery in action, we study 
a high-significance anomaly in the WMAP data, which 
resembles the effects of a modulation of the primordial 
power spectrum. Explicitly, one can construct estimators 
based on the covariance for a fc-space modulation of the 
primordial power spectrum V x (k) with the form 

V x (k)=V x (k)[l+g(k)}. (23) 

If we take the bi-Copernican hypothesis that the universe 
has no preferred orientation then the expectation is that 
g(k) = 0. However current analyses of the WMAP maps 
strongly favor a model in which <?(k) has quadrupolar 
components 

5 (k) = J2 92MoY 2M (k). (24) 

\M\<2 



Furthermore, the preferred g 2 M are planar, with g 2 M oc 
<5mo in ecliptic coordinates [HI [T7J [19] : 

.9(k)=ff 2 o^yj(3cos 2 ^-l), (25) 

where 9^ is the angle from the ecliptic pole. This form 
resembles the model proposed by [20] (and other authors, 
e.g. [21]) which was the motivation for original detection 
made in [25]. A missing factor in the original version of 
[20] resulted in the ecliptic orientation being obscured; 
however, this was corrected by [15j . The ecliptic align- 
ment of the detected effect strongly suggests an instru- 
mental systematic or solar-system origin. The signal is 
present at 9<r in the W band [TO], but varies strongly 
between detectors at the same frequency [15] [17] , which 
singles out an instrumental explanation, although [19] 
have also checked the contribution of zodiacal light and 
found a negligible effect. 

Here we will continue the work of [15] , using optimal 
quadratic maximum-likelihood estimators to study the 
primordial modulation effect. These estimators are often 
favorable to the Gibbs-sampling approach of [THl [22] for 
their speed, and the ease with which they can be modified 
to test various systematic effects [HI [17] . In the current 
application, the quadratic estimator compresses millions 
of correlations between thousands of observed modes to a 
small handful of parameters, and should be effectively in- 
distinguishable from an exact likelihood analysis [17) . In 
the formalism of the previous section, the quadratic es- 
timator for <?2M has S = sj = 0, and the weight function 
can be written as [HI [23] 

W(h,l 2 ,2) = T_ C hl2 , (26) 

where the Ci x i 2 matrix is given by 

C hh = 4tt / dlnkV x (k)A h (k)A l2 (k). (27) 

The A;(fc) used here are the angular CMB transfer func- 
tions. 

The largest expected instrumental effects which can 
produce ecliptic-aligned anomalies are inhomogeneous 
pixel noise levels and beam asymmetries. The pixel- 
uncorrelated component of the instrumental noise is al- 
ready accounted for in current analyses, and it has been 
argued that these estimators are insensitive to percent- 
level changes in the noise amplitude [12]. We agree with 
this result: the mean-field for the quadratic estimator of 
primordial power modulation due to WMAP noise inho- 
mogeneities is less than la for all V- and VF-band dif- 
ferencing assemblies (DAs), and so percent-level changes 
in the noise level do not have appreciable effects. In the 
bipolar power spectrum formalism, this is because inho- 
mogeneous instrumental noise produces coefficients of the 
form Af^ oc const, (in the notation of Q~7]), which more 
closely resembles a modulation of the observed power 
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spectrum in real space than a modulation of the pri- 
mordial power spectrum in k space [17j . Anomalies are 
also seen in estimates of g 2 M formed from cross corre- 
lations between maps with different noise realizations, 
which suggests that noise cannot be the dominant effect 

mm- 

It has been argued that beams must provide at least 
a significant source of bias for estimates of g(k) [15], if 
not a complete explanation [17], although [19] have also 
studied the effects of beams and concluded that they are 
unimportant. Here we will address this issue using the 
new tools of the previous section. We obtain the coeffi- 
cients of the instrumental beam directly from the WMAP 
five-year published beam maps |24j by a brute-force dis- 
crete harmonic transform, with beam center determined 
simply from the maximum pixel. For each beam map, 
we then scale the resulting B lm such that B 00 = 1, and 
average the A- and -B-side beams together, as appropri- 
ate for simulating the effect of differencing on the final 
map [H] . Finally, we scale these averaged beam transfer 
functions at each I so that the Bi components agree with 
the published WMAP transfer functions. These are de- 
rived from the same TOD observations of Jupiter which 
are used to create the beam maps, but do not suffer from 
pixelization effects. Without this scaling, our Bi would 
still agree with the published values at a level of better 
than 1% for I < 600. 

To get a feeling for the expected effects we begin by 
evaluating Eq. (22). We only calculate the bias due to 
the A^ 1 ) terms of the beam covariance, however we will 
later verify numerically that the A^ 2 ^ terms are not sig- 
nificant for this application, as expected. As we are ul- 
timately only interested in the low-Z multipoles of the 
beam mean-field it is sufficient to have a model for the 
scan strategy swlm at low-/, and for correspondingly 
small 5*. We will initially use s ma x = 6. To calculate the 
swlm we use the analytical method of [7 , outlined in 
Appendix [X] which provides an excellent approximation 
to the true scan strategy on such large scales. Following 
the notation there we use a toy model for WMAP, with 
a spin period of two minutes, a precession period of one 
hour, and scan angles of 9b = 70°, 9 p = 22.5°. These are 
design values for the fiducial center of the WMAP focal 
plane, and were achieved with good accuracy in flight 
[2"5] . The position of the individual detectors within the 
focal plane does effectively vary 9b- We have not cor- 
rected for this; however, as we will discuss, this does not 
significantly affect our results. The expected mean-field 
biases are presented in Table [I] for the V- and VF-band 
DAs. We initially limit our analysis to Z max = 400 for 
comparison with earlier results. The inverse-variance fil- 
ters use isotropic noise with the appropriate power spec- 
trum. This evaluation takes only a few seconds for such 
low s max . The fast precession of the WMAP spin axis 
gives the scan-strategy azimuthal symmetry in ecliptic 
coordinates, which makes swlm <x Saw- Because the 
bias is proportional to this quantity, it also has this az- 
imuthal symmetry, which would explain the planar struc- 



ture of the detected modulation, and the alignment of the 
detected effects with the ecliptic poles. The north-south 
symmetry of the scan strategy also restricts s w lm to 
even-L. The mean-fields for I = 2 are predicted to be 
large, and detectable at many sigma. Higher multipoles 
(I = 4, 6) receive much less significant contributions from 
beams, and are also not observed to be anomalous in the 
data [15] . Thus, beams seem to be a likely explanation 
for the detected anomaly. 

Before we move on to the analysis of the WMAP data 
itself, we consider some of the insights which the analytic 
approach makes possible. Because the detected anomaly 
is quadrupolar, it depends on the scan strategy only up 
to S = 2 (as there are no L = 2 modes for larger S). 
This means that if beams are the explanation they can 
only be sourced by the beam's dipole and ellipticity com- 
ponents. These modes are well constrained by the beam 
maps, so we expect our calculations to be quite accu- 
rate. As already mentioned, the effective 9b which we 
use for our scan-strategy calculation differs for some de- 
tectors due to their position in the focal plane. We find 
that d(2W2o)/d9b = 2%/deg. about the fiducial value. As 
all of the V- and PF-band detectors are clustered within 
1° in azimuth of the center of the focal plane, we ex- 
pect at most 2% errors in the scan-strategy coefficients. 
This corresponds to < 0.2er effects for the biases which 
we have derived. The analytical calculations also reveal 
that the dominant contribution to the quadrupolar bias 
is given by the S = 2 modes, rather than 5 = 1. For 
our purposes, this is quite fortuitous- there are a num- 
ber of opportunities in this calculation for 180° errors in 
the beam orientation angle, however any such errors will 
not significantly affect our results. 

Another question which can be asked is whether it 
might be possible to verify the effects of beams with 
an estimator more optimally designed to detect them. 
A Fisher-matrix calculation shows that the g 2 Ai estima- 
tor has a typical correlation coefficient of 0.9 with the 
optimal quadrupole "scan-strategy" estimator with the 
known WMAP beams, and so the g 2 M estimator is effec- 
tively optimized to detect beam effects. 

Finally, it is interesting to note that the bias for the Q- 
band data has the opposite sign of that for the V and W 
bands. This effect has already been observed in the data 
by [19], although without explanation. It is due to the 
fact that the semimajor axis of the beam ellipticity for the 
Q band is oriented parallel to the scan direction, while 
the V/W axes are perpendicular [25]. A 90° rotation 
corresponds to a sign flip for the S = 2 modes of the 
beam which dominate the bias. This feature provides 
strong evidence that the dominant effect which sources 
the quadrupolar effect is beam asymmetries. 

We now turn to the question of whether beam ef- 
fects are sufficient to explain completely the quadrupole 
anomaly, as seems likely from the significance levels in 
Table [I] It is straightforward to fold asymmetric beam 
effects directly into our analysis [TH] of the WMAP five- 
year data [21]; we have not yet upgraded to the seven- 
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DA ( ffa0 ) (cr) (<?4o) {&) (flreo) (<r) 



Ql 

Q2 

VI 

V2 

Wl 

W2 

W3 

W4 



-0.33 (12.1) 
•0.33 (12.3) 
0.17 (6.51) 
0.17 (6.74) 
0.27 (9.10) 
0.31 (9.79) 
0.33 (9.99) 
0.27 (8.63) 



0.030 (0.83) 
0.029 (0.81) 
0.031 (0.86) 
0.032 (0.92) 
0.043 (1.07) 
0.042 (0.97) 
0.037 (0.85) 
0.045 (0.95) 



-0.003 (0.07) 
-0.003 (0.08) 
-0.003 (0.07) 
-0.002 (0.06) 
-0.002 (0.05) 
-0.003 (0.06) 
-0.002 (0.05) 
-0.003 (0.05) 



TABLE I: Analytic predictions for the beam mean-field bias 
to the primordial power asymmetry estimator, with Z max = 
400. The significance (a) is given by the mean-field divided 
by the estimator noise. 



year data, which are very consistent |17j . In the quadratic 
estimator formalism which we use here, the estimator 
mean-field and Fisher matrix are determined on the cut 
sky by Monte Carlo. By incorporating beam effects into 
the CMB simulations using Eq. the Monte-Carlo 
mean-field will include the contribution from beams. 
Note that this approach will remove the subdominant 
contribution due to the A' 2 ) terms, as well as the A^ 1 ) 
terms which we used for the analytical calculation, al- 
though we have verified numerically that the A^ 2 ) terms 
do not contribute more than a few percent to the beam 
mean-field. For the convolution we use s max = 2, as higher 
terms do not contribute to the quadrupolar anisotropy. 
To obtain a minimum- variance estimator, we should also 
incorporate beam effects into the inverse-variance filter- 
ing operation, as discussed in the previous section, how- 
ever we find that the estimator noise variance for the 
asymmetrically convolved simulations is less than 10% 
higher than that without beam effects, and so this im- 
provement would not have a significant effect on our re- 
sults. 

The significance of the measured <?2M is plotted in 
Fig. [T] for the foreground-reduced WMAP V- and W- 
band DAs, incorporating all of the usable signal in the 
WMAP data by taking Z max = 1000. It can be seen that 
the mean-field subtraction of beam effects results in data 
which are consistent with the isotropic model. The non- 
Smo values are very similar for all detectors, as one would 
expect for an isotropic sky, given that a large fraction of 
the estimator "noise" is due to the CMB fluctuations 
themselves, which are common between detectors. The 
variation for the Smo modes is more significant, indicat- 
ing that there may be some small errors in the mean-field 
subtraction. The measured value of g2o is not significant 
in any detector, however, with the exception of W4 which 
shows a large negative bias even after subtraction of the 
beam mean-field. It should be kept in mind that with- 
out beam subtraction, each DA shows 6 — 9er effects in 
the Smo modes. We believe that the residual anomaly in 
W4 may be attributed to the effects of correlated noise. 



The W4 DA has clearly correlated noise even after the 
prewhitening stage of the WMAP analysis [25] . Its 1// 
knee frequency is several times larger than any of the 
other DAs [57]. A negative bias to c/20 from correlated 
noise is also expected analytically. For subdominant cor- 
related noise, the effects of striping can be modeled as a 
convolution of the pixel-uncorrelated instrumental noise 
with a narrow beam which has its semimajor axis along 
the scan direction. As we have already seen with the Q- 
band data, this leads to a negative bias in (720 • We can 
further test this hypothesis by forming our quadratic es- 
timator from pairs of maps with uncorrelated noise. For 
this we use the Wi data for individual years. Autocor- 
relating the data for any single year and correcting for 
beam effects, we continue to find a large negative bias in 
520 (albeit with slightly less statistical significance due 
to the larger instrumental noise for only a single year of 
data). Cross-correlating data from any pair of separate 
years, however, the effect disappears. Although we do 
not plot them here, the Q-band data are also completely 
consistent with 520 = after correction for beam effects. 
We therefore assert that beam effects provide a sufficient 
explanation for the detected anomaly. Based on the av- 
erage diagonal elements of our simulation Fisher matri- 
ces, we place a conservative limit on any single mode of 
|<?2mI < 0-07 at 95% confidence. The corresponding limit 
for the power spectrum of g^M is Cf 9 < 0.003. 

Disagreement with other results. - In [15j we tested the 
effect of beam asymmetries using the simulations of pQ , 
which appeared only partially to explain the strong de- 
tected signal. However, on closer inspection, we find 
problems with these simulations, indicating that they are 
not representative of the effect of the WMAP beams. 
Applying the modulation estimator to these simulations 
with no noise, we find that they show a strong mean- 
field, with a spatial pattern which is identical between 
detectors, as expected analytically. The mean-field does 
not, however, have the purely planar structure associated 
with a Smo pattern in ecliptic coordinates. Such an er- 
ror is most likely due to errors in the beam orientation 
angles a:,; however, private communication with Eriksen 
et. al. has been unsuccessful in revealing the precise ori- 
gin of this discrepancy. Our results are also discrepant 
with |19) . who analyzed one of the simulations from [1] 
and found no beam effects. We believe this is due to two 
factors: (i) the simulation which they analyzed is one in 
which we also saw small effects in our original work [15] ; 
and (ii) the nonplanar structure of the mean-field in the 
simulations poses difficulty for the estimator used by [19j , 
which searches explicitly for the planar mode. 



V. CONCLUSIONS 

Beam asymmetry effects can be more important for 
anisotropy estimators than for power spectrum analysis, 
because their effects generally enter at lower order in the 
asymmetry. Beam asymmetries fit nicely into the larger 
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FIG. 1: Significance of the observed WMAP primordial- 
power-moduiation quadrupole, with correction for beams, 
for the WMAP V-band foreground-reduced data, limited to 
imax = 1000. This is essentially a beam-corrected version of 
Fig. 9 in [15]. The gray histograms are from the 500 simula- 
tions which are used to determine the estimator Fisher matrix 
and mean-field. The dashed vertical lines are for the observed 
data. Detailed interpretation is provided in the text. 



formalism of quadratic anisotropy estimators. They re- 
sult in a mean-field bias which directly traces the scan 
strategy s w lm, and can be calculated analytically on the 
full sky or determined from Monte-Carlo simulations on 
the cut sky. 

Beam effects appear to provide a sufficient explanation 
for the 9cr detection of an apparent quadrupolar modu- 
lation of the primordial power spectrum in the WMAP 
data. We note that the WMAP team already incorpo- 
rates the effects of beam asymmetry into their power 
spectrum analysis (where it is a much smaller effect in 
any case), and so the resolution of this anomaly should 
not have any effect on their cosmological parameter con- 
straints. 

All of this work will apply directly to the Planck exper- 
iment, which has a less symmetrizing scan strategy than 
WMAP. Planck's increased sensitivity also opens up the 
field for the precision analysis of interesting astrophysical 
secondaries such as the anisotropic signal from gravita- 
tional lensing. The tools and techniques which we have 
discussed here may also be extended straightforwardly 
for use with polarization data. 
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Appendix A: Toy Scan Strategy 

Here we consider a simple model of a scan strategy 
which serves as a good approximation to typical satellite 
experiments. It consists of: 

(1) a beam at an angle 9^ to the satellite spin axis, which 
rotates with period t s ; 

(2) a precession at an angle 9 P to the antisolar direction, 
with period r p ; and 

(3) a continuous repointing of the antisolar direction as 
the observer orbits the sun. 

If t s <C t p <C lyear, then w(Cl p ,s) can be calcu- 
lated analytically [7J. First, we calculate the quantity 
v(Cl p ,s) = J2i£ P elsai an d th cn we form w(Cl p ,s) — 

v(n p ,s)/v(n p ,o). 

To calculate v(Q p ,s), begin in a coordinate system 
which places the spacecraft spin axis along the +z axis. 
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Rotation about the spin axis in these coordinates gives 
vW(f2 p , s) oc 5(6 — 9b)e* s0 . Expanding this using the ap- 
propriate spin harmonics gives 



[v^(n p , s )] lm = Ks m0s Y lo (e b ,o), 



(Al) 



where K is some constant. We can then rotate out to 
place the precession axis along the +z axis, obtaining 
v(p, s) in precession coordinates: 

{v ( v (n p , s )] lm = ^ D l mm „ (o, e p , o) [v<u (n p , 8 )\ lm „ . 

to" 

(A2) 

Rotation about the precession axis removes all but the 
m = components giving 

[vW(n p ,s)] lm = 5 m0 D l 00 (0,d p ,0)[v^(Q p ,s)] w 

= 5 mQ KPi(cos9 p ) s Y l0 (9 b ,0). (A3) 

Finally, the precession axis is rotated 90° down to the 
ecliptic plane, and again only m = modes in the new 
coordinates are taken, to effect the azimuthal averaging 
given by the yearly rotation about the sun: 

[v(Qp, s)] lm = 5 m0 KPi(0)Pi(cos9 p ) s Yi (9 b , 0). (A4) 

This can be used to calculate w(Q p , s) in ecliptic coordi- 
nates. Note that only multipoles with I even are nonzero 
due to the north-south symmetry of the scan pattern. 
This is enforced by the dependence on P/(0). 



Appendix B: Beam asymmetries and power spectra 

For the CMB temperature power spectrum, beam 
asymmetries are only important at high I (i.e. below the 
beam scale). On such scales, a pseudo-C; analysis is usu- 
ally employed in which a weight function is applied to 
the observed sky and the empirical (pseudo) power spec- 
trum of the weighted sky is taken. The expectation value 
of the pseudo-C; , after removal of the bias due to instru- 
ment noise, is linearly related to the true power spectrum, 
Cf e . In this Appendix, we calculate this relation in the 
presence of beam asymmetries; for related work see [IT] . 
Good performance can be obtained from pseudo-C; esti- 
mators with a careful choice of weight function [2§] , such 
as a local approximation to the optimal inverse signal- 
plus-noise weighting. Below the beam scale, the signal 
is exponentially suppressed and weighting by the inverse 
variance of the pixel noise is close to optimal. 



We noted in Sec. [TT] that, for full-sky observations and 
uniform weighting of the data, beam asymmetries only 
affect the power spectrum at second order. This is no 
longer true with anisotropic weighting, which generally 
will arise from the survey geometry or inhomogencitics 
in the noise. In this case, we generalize Eq. Q to in- 
clude a spin-0 pixel weight function on the right-hand 
side. Equation (|8| then still holds, but oWlm is no longer 
necessarily zero for L ^ 0. Writing the integral of three 
spin harmonics in terms of 3j symbols, we have 



LAIS Im 

( V I l \ ( V I L 
< \ -m' m M I l -S S 



(2/ + l)(2Z' + l)(2L + l) 



47T 



BlS-SWL M @lm- ( B1 ) 



Inserting this in the definition of the pseudo power spec- 
trum, 



Ci 



21 



— y |e 
+ 1 z^i 



.2 

Ira | j 



(B2) 



and taking the expectation value, we find 



- ^ (2/ + l)(2£ + l) (V I L\(V I L\ 

{ v) = ks ^ {o-ss){o-s> S'j 

xB^B^^s.s^WLCf , (B3) 

where we have used an orthogonality relation for the ij 
symbols. Here, the scan strategy and weight function are 
encoded in the cross spectra 



(ss')Wl 



2L 



— E 

+ 1 V 



swlms'Wlm- 



(B4) 



For mildly asymmetric beams, or a wide distribution of 



crossing angles, the sums over S and S' in Eq. (B3| can 
be truncated after only a few terms. In this case, the 
effect of beam asymmetries on the mean pseudo-C; can 
be computed efficiently with no further simplifying as- 
sumptions. 

For the case of symmetric beams (Bis = $SoBio), or 
for a uniform distribution of observation angles in each 
pixel (swlm = s w lmSsq), Eq. (B3) reduces to the usual 
result [29] for symmetric beams. 
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